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Abstract 

Quark-gluon decays of heavy paraquarkonia 1 Sq{QQ ) —> qqg are 
investigated with account of the masses of final quarks. The decay 
widths and the energy distributions of the final quarks and gluons are 
calculated in dependence on the relative quark masses. The strong 
collinear enhancement of the gluon energy distribution at the end of 
the spectrum is shown to take place in all such decays of r/ c and '% 
mesons except the decay rjb —> ccg. The total decay width is shown 
to have an essential dependence on the final quark masses. The cor¬ 
responding branching ratios of g c and rjb mesons are numerically es¬ 
timated with a good agreement of Br(r] c —► ssg) with experimental 
data on r/ c decays. 
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Many-partical decays of heavy quarkonia such as 1 So(QQ ) —> 3g,ggg, 
3 Si(QQ) —> 4g,gggg depend on the 3 g- and ggg-inter act ions of the final 
quarks and gluons and can give the useful information about these interac¬ 
tions. In particular, these interactions give rise to the specific thrust |], @ 
and invariant masses |3], Q distributions, to the acomplanarity of the four- 
particle decays ||, to the collinearization effect of final gluons 0, etc. It 
should be noted, however, that such effects can be essentially affected by 
masses of final quarks. For example, the four-particle quark-gluon decays of 
the orthobottomonium with the production of uu-, dd- or ss-pairs exhibit a 
collinear enhancement whereas this effect in such decays with the production 
of cc-pair is absent because of the relatively large mass of c-quark |7j . Hence, 
the effect of the final quark masses can be essential in the decays of heavy 
quarkonia with production of gg-pairs and it should be taken into account in 
investigations of such decays. 

In this letter the quark-gluon decay of heavy paraquarkonium 1 S 0 (QQ ) —> 
qqg is investigated taking account of the masses of final quarks. The differ¬ 
ential and total decay widths are calculated in tree approximation in depen¬ 
dence on the final quark masses and are discussed in comparison with the 
experimental data on p c decays. The branching ratios of rjb decay are also 
predicted and discussed. 

The decay of a heavy paraquarkonium into a quark-antiquark pair and 
a gluon 1 S 0 (QQ) —>► qqg is described by two graphs shown in Fig. |]. The 
amplitude of this process in the limit of static quarks in the quarkonium can 
be presented in the form: 
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(ta)ap 9s^( 0 ) UagF a V) 
2VN C 2m^/uj£ l £ 2 p 2 (Vk) 


( 1 ) 


where t a are the generators of the color group SU(N C ), a = 1,2,, AT? — 1 
and a, (3 = 1,2 ,..., N c are the color indices, g s is the strong charge, i/j( r) 
is the nonrelativistic wave function of a paraquarkonium in the coordinate 
space, lu, £i and e 2 are the energies of the gluon, quark and antiquark re¬ 
spectively, (j^aS = {u a (p 1 )'y /1 u /3 (-p 2 )) is the final quark current, = 
e^i /p ak p e^, e“ and k M are the polarization and wave vectors of the gluon, 
p = pi + p 2 , pin and p 2fl are four-momenta of the final quark and antiquark, 
Vf, — (2m, 0) is the four-momentum of a paraquarkonium in its rest frame, 
and m is the mass of the heavy quark in the quarkonium. 
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The differential probability of the decay after summation over colors and 
polarizations of the final particles can be written in the form: 
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where F c = (N 2 — 1)/(8 N c ) is the color factor of the SU(N C ) group (in QCD, 
N c = 3 and F c = 1/3), a s = g 2 /(4n) is the strong coupling constant, and 
/i = m q /m is the relative mass of a final quark (antiquark). 

The probability (^) of the three-particle decay 1 S 0 (QQ) —>► qqg of the 
paraquarkonium depends on two independent variables and as a function of 
two independent energies can be presented in the form: 
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where y\ = £\/m, y 2 = £ 2 /m and x = u/m are the relative energies of the 
quark, antiquark and gluon respectively satisfying the energy conservation 
low x + yi + y 2 = 2. The expression (j3[) is rather simple and allows for 
the relative mass of the final quark. In the particular case of g = 0 it is 
consistent with the corresponding result presented in the more complicated 
form in Ref. |[|| . 

Integrating Eq. (E[) over the momenta of the quark pi and the antiquark 
P 2 or over the momenta of the antiquark P 2 and the gluon k, we obtain the 
distribution in the energy x of the final gluon or that in the energy y\ of the 
final quark in the form: 
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For further analysis it is convenient to normalize these distribution func¬ 
tions as 
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where 


r 2s = rCSc(QQ) - 2 g) = F c ^ |V>(0)| 2 

TTT 


is the two-gluon decay width of the paraquarkonium. The normalized dis¬ 
tribution functions f g and f q in the gluon and quark energies are presented 
in Figs. [| and [| at the relative final quark masses // = 0.33, // = 0.10 and 
H = 0.03 corresponding to the decays —> ccg, r/ c —> ssg and rjb —> ssg 
respectively. 

As seen from Fig. ^ the gluon energy distribution function f g has the 
strong enhancement at the end of the spectrum in the case of the small 
relative mass // = 0.1 {jj c —> ssg) (as well as in the case of n — 0.03 ( 7j b —> 
ssg)), whereas this enhancement is absent in the rjb —> ccg decay because of 
the large relative mass of c-quark (/i = 0.33). The origin of this effect is the 
quark-antiquark collinear enhancement. 

The quark energy distribution function f q (see Fig. |3|) has a small en¬ 
hancement for /i = 0.10 and /j = 0.03 caused by an emission of the soft 
gluons. 

Integrating the distribution (|4j) over the gluon energy x (or the distribu¬ 
tion (H) over the quark energy yi) we have obtained the total width Y qqg of 
the quark-gluon decay of the paraquarkonium with account of the masses of 
final quarks in the form: 


- | (4 - !?) AT?} . (6) 

The dependence of this width on the masses of the final quarks can be de¬ 
scribed by the function 
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shown in Fig. [|. The points on the curve in Fig. [| correspond to the rj c meson 
decay rj c —> ssg and to the rj b meson decays rj b —> ccg and rj b —> ssg. As seen, 
the dependence of the width T qqg on the hnal quark masses is essential in 
analyzing the decays of the rj c and rjb mesons. 

For numerical estimations we use further the values a s (rrib) ~ 0.16 and 
a s {m c ) ~ 0.21 resulting from a s (mz) = 0.118±0.003 0. Using the width (|^) 
and the theoretical expression for the hadronic width of the paraquarkonium 
Thad = r 2g (l — Cds/it ) with the C = 5.84 for charmonia and C = 5.41 for 
bottomonia (l|] we have calculated the branching ratios Br = T qqg /T j ia d of 
1 So(QQ) -^ QQ9 decays. The results are presented in Table |T| for the hnal 
quark masses varying from their current values (the first values of m q and 
Br) to constituent ones (the last values of m q and Br). 

The result concerning the rj c —> ssg decay can be compared with the 
experimental value of Br(g c —> KK+X) (where X is anything else) presented 
in Table |]. The theoretical prediction of Br(// C —> ssg ) = 12% at the current 
hnal quark masses is in good agreement with the experimental value Br(r/ C —> 
KK + X) = 11.6 ± 2.4%, i.e. the summed width of the rj c —> KK + A" type 
decays is saturated mainly by the subprocess g c —> ssg. One may hope 
that the branching ratios Br(r/;, —> ccg ) ~ 2.5% and Br(r/ & —> ssg ) ~ 13% 
presented in Table |] also give a good approximation for the branching ratios 
of the parabottomonium decays )]b —> DD + X and rj b —> KK + X. 

We summarize the results of this letter as follows: 

The quark-gluon decays of heavy paraquarkonia 1 S 0 (QQ) —> qqg are in¬ 
vestigated in the tree approximation taking account of the hnal quark masses. 

The distribution functions in quark and gluon energies are obtained and 
discussed. The distribution function in the gluon energy is shown to have 
the strong collinear enhancement in all the decays of rj c and rjb mesons of 
1 Sq(QQ) —> qqg type except the decay rjb —■► ccg. 

The total width of 1 Sq{QQ ) —> qqg decay is calculated and its dependence 
on the hnal quark masses is shown to be essential for decays of rj c and rjb 
mesons. The corresponding branching ratios of rj c and rjb mesons are numer¬ 
ically estimated with a good agreement of the branching ratio Br(?; c —> ssg ) 
with experimental value of Br (rj c —> KK + A"). 
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Table 1: Theoretical predictions for the branching ratios Br = T^/Thad of 
the quark-gluon decays of heavy paraquarkonia. 


Decay 
Vc ssg 
Vc —> ddg 
Vc uug 
Vb -> ccg 
r] b -> ss# 
rj b —► ddg 
lib —> uug 


m q (MeV) 
150 4 - 450 
7 4 - 300 
4 4- 300 
1500 

150 4- 450 
7 4- 300 
4 -r 300 


Br (%) 
12^5 
35^7 
38 4- 7 
2.5 

13-5-8 
27 -5- 10 
30 - 5 - 10 
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Tabic 2: Experimental values || of the fractions of the r] c 
decays (.X is anything else). 


Modes 

Fractions 

KKn 

K*°K~tt + + c.c. 
K + K-Tr + n- 
2 (K+K~) 

(5.5 ±1.7)% 
(2.0 ± 0.7)% 

(2.o±8:S)% 
(2.1 ±1.2)% 

KK + X 

(11.6 ±2.4)% 


KK ± X type 
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Figure captions 


Fig. 1. Graphs describing the quark-gluon decay of a heavy paraquarkoninm 

1 S 0 (QQ) -> qqg- 

Fig. 2. Distribution functions f g in the gluon energy x at the relative masses 
of final quarks g = 0.33 (rg —> ccg) and g = 0.10 (r/ c —> ssr/). 

Fig. 3. Distribution functions f q in the quark energy y\ at the relative masses 
of final quarks g = 0.33 (■% —> ccg), g = 0.10 (r/ c —> and g = 0.03 

(Vb ssg). 

Fig. 4. The normalized width / tot of the decay 1 S 0 (QQ) —> qqg as a function 
of the relative final quark mass g. 
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